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Abstract. Eternal solutions of parabolic equations (those which are defined 
for all time) are typically rather rare. For example, the heat equation has 
exactly one eternal solution - the trivial solution. While solutions to the heat 
equation exist for all forward time, they cannot be extended backwards in 
time. Nonlinearities exasperate the situation somewhat, in that solutions may 
form singularities in both backward and forward time. However, semilinear 
parabolic equations can also support nontrivial eternal solutions. This article 
shows how nontrivial eternal solutions can be constructed for a semilinear 
equation that has at least two distinct equilibrium solutions. The resulting 
eternal solution is a heteroclinic orbit which connects the two given equilibria. 



1. Introduction 

Consider 



(i) ^ = M(, I ) + ^ ttl (xK(t,x), 



at 

where ( £ 1 and x € K n , and the eij are bounded and smooth. In this article, we 
consider eternal solutions, those classical solutions u which satisfy ([T]) for all time 

t e m. 

This kind of equation provides a simple model for a number of physical phenom- 
ena. First, choosing the right side to be Au — u 2 + a\u results in an equation which 
can represent a model of the population of a single species with diffusion and a 
spatially- varying carrying capacity, a±(x). As a second application, this equation is 
a very simple model of combustion. If a\ is a positive constant, then the equation 
supports traveling waves. Such traveling waves can model the propagation of a 
flame through a fuel source. 

Equations of the form ([1]) have been of interest to researchers for quite some 
time. Existence and uniqueness of solutions on short time intervals (on strips 
(—to, to) x R") can been shown using semigroup methods and are entirely standard 
[T5] , However, there are obstructions to the existence of eternal solutions. Aside 
from the typical loss of regularity due to solving the backwards heat equation, there 
is also a blow-up phenomenon which can spoil existence in the forward-time solution 
to ([I} . Blow-up phenonmena in the forward time Cauchy problem (where one does 
not consider t < 0) have been studied by a number of authors [Hj [5] [13] [H] [2] [T7] 
[TS] , More recently, Zhang et al. ([H] [H] [T3]) studied global existence for the 
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forward Cauchy problem for 




.p 



V{x)u 



for positive u, V. Du and Ma studied a related problem in [4] under more restricted 
conditions on the coefficients but they obtained stronger existence results. In fact, 
they found that all of the solutions which were defined for all i > tended to 
equilibrium solutions. 

Eternal solutions to (p} are rather rare. Most works which describe blow-up 
make the assumption that the solution is positive. Unfortunately, blow-up is much 
more difficult to characterize in the general situation, and understanding exactly 
what kind of initial conditions are responsible for blow-up in the Cauchy problem 
for |T|) is an important part of the question. 

As an aside, the boundary value problem that results from taking x £ 57 C K ra 
for some bounded f2 (instead of x € has also been discussed extensively in the 
literature [7] [8] [3] . For the boundary value problem, all bounded forward Cauchy 
problem solutions tend to limits as \t\ — > oo, and these limits are equilibrium 



The existence of eternal solutions is a difficult problem, because the backward- 
time Cauchy problem is well known to be ill-posed. Obviously, equilibrium solutions 
are trivial examples of such eternal solutions, and in [lOj it was shown that they can 
exist. It is not at all clear that there are other eternal solutions, and indeed there 
may not be. In this article we assume the existence of a pair of nonintersecting 
equilibrium solutions and construct a heteroclinic orbit which connects them. (A 
heteroclinic orbit is a special kind of eternal solution, whose limits as t — > ±oo are 
equilibria.) 

For simplicity and concreteness, we will work with the more limited equation 



where is a smooth function which decays to zero. We follow the general tech- 
nique for constructing "ancient solutions," which was used in a different context 
by Hamilton and Perelman. It should be emphasized that the technique exam- 
ined in this article can handle the problem for (|TJ) in full generality (under mild 
decay assumptions for the coefficients aj), though this complicates the exposition 
needlessly. 

This simpler model still provides insight into applications, as it is still a model 
of the population of a single species, with a spatially- varying carrying capacity, <j>. 
Indeed, one easily finds that under certain conditions the behavior of solutions to 
is reminiscent of the growth and (admittedly tenuous) control of invasive species 
P]. It is the control of invasive species that is of most interest, and it is also what 
the structure of the space of heteroclines describes. In one of the examples given 
in [TU], there is one "more stable" equilibrium, and several other "less stable" ones. 
The more stable equilibrium can be thought of as the situation where an invasive 
species dominates. The task, then, is to try to perturb the system so that it no 
longer is attracted to that equilibrium. An optimal control approach is to perturb 
the system so that it barely crosses the boundary of the stable manifold of the the 
undesired equilibrium, and thereby the invasive species is eventually brought under 
control with minimal disturbance to the rest of the environment. Such an optimal 
control approach, though, is beyond the scope of this article. 



solutions. 



(2) 



du(t,x) _ d 2 

at ~ 



: ' X) 



dx 2 



u 2 (t, x) + 4>(x), 
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2. EQUILIBRIUM SOLUTIONS 

We choose </>(x) = (x 2 — 0.4)e~ x2 / 2 . It has been shown in [10] , that in this 
situation, there exists a pair of equilibrium solutions / + , /_ with the following 
properties: 

(1) /+ and /_ are smooth and bounded, 

(2) /_)_ and /_ have bounded first and second derivatives, 

(3) /_)_ and /_ are asymptotic to 6/x 2 for large x, and so both belong to L 1 (IR), 

(4) f+{x) > f_(x) for all x, 

(5) there is no equilibrium solution f% with f+(x) > f%(x) > f-{x) for all x, 
and additionally, there exists a one-parameter family g c of solutions to 

(3) Q = g'^x)-gl{x)+<j> c {x) 
with 

(1) ce [0,1), 

(2) 5o = /- and O = 4>, 

(3) 4> a {x) < 4>b(x) and g a (x) > gb(x) for all x if a > b. 

The latter set of properties can occur as a consequence of the specific structure 
of For instance, consider the following result. 

Proposition 1. Suppose /_ G C 2 ' Q (R) satisfies the above conditions and addition- 
ally, there is a compact K C K with nonempty interior such that /_ is negative on 
the interior of K and is nonnegative on the complement of K. Then such a family 
g c above exists. 

Proof. (Sketch) Work in T/_C 2 ' Q (IR), the tangent space at /_. Then © becomes 
its linearization (for h c , say), namely 

(4) = t%(x)-2f-(x)h c (x) + (<l> c -<j>). 
Consider the slighly different problem, 

(5) - y"(x) - 2f_(x)y(x) + v(x)y(x), 

where v is a smooth function to be determined. If we can find a v < such that 
y > and y — ^ as — !• oo, then we are done, because we simply let vy = 4> c — 4> 
in (f?]). In that case, h c = y has the required properties. We sketch why such a v 
exists: 

• If v = 0, then y = is a solution, giving g c = /_ as a base case. 

• If v(x) = — 2||u|| 00 /3(x) for (3 is a smooth bump function with compact 
support and (3\K — 1, then the Sturm-Liouville comparison theorem implies 
that y has no sign changes. We can take y strictly positive. However, in this 
case, the Sturm-Liouville theorem imples that there are no critical points 
of y either, so y may not tend to zero as |x| — * oo. 

• Hence there should exist an s with < s < 2||it||oo such that if v(x) = 
— s(3{x), then y has no sign changes, one critical point, and tends to zero as 
|a;| — > oo. This choice of v is what is required. (The precise details of this 
argument fall under standard Sturm-Liouville theory, which are omitted 
here.) 

□ 
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In what follows, we shall not be concerned with the exact form of </>, but rather 
we shall assume that the above properties of the equilibria hold. Many other choices 
of (f> will allow a similar construction. 

Lemma 2. The set 

(6) W = {v£ C 2 (R)\f-(x) < v(x) < f+(x) for all x} 

is a forward invariant set for @. That is, if u is a solution to @ and u(fo) € W , 
then u(t) G W for all t > t . 

Proof. We show that the flow of @ is inward whenever a timeslice is tangent to 
either /_ or / + . To this end, define the set B 

B = {d£ C 2 (R)\f^(x) < v(x) < f+(x) for all x, and there exists an x 

such that v(xo) = f+(x) or v(xq) = /_(a;)}. 

Without loss of generality, consider a v E B with a single point of tangency, v (xo) — 
f-(xo). At such a point xo, the smoothness of v and /_ implies that Av(xq) > 
Af-(xo) using the maximum principle. Then, if u is a solution to (U) with u(0, x) — 
v(x), we have that 

= Av(x )~v 2 (x ) + ^x ) 



hence the flow is inward. One can repeat the above argument for each point of 
tangency, and for tangency with / + as well. □ 



Lemma 3. Solutions to the Cauchy problem 

\u(o,x) = u(x) e w c 



(7) 
where 

W c = {ve C 2 {R)\g c {x) < v{x) < f+{x) for all x} 

for c £ [0, 1) have the property that they lie in L 1 n for all t > 0. We shall 

assume that U has bounded first and second derivatives. 

Additionally, when c e (0, 1), solutions to ([7]) cannot have /_ as a limit as 
t —>■ oo. 

Proof. The fact that solutions lie in L 1 nL°° (R) is immediate from Lemma[2]and the 
asymptotic behavior of /+, /_ (Section 4 of [TO])- Observe that for each c € [0, 1), 
W c is forward invariant, and that W a C W\, if a > b. Since /_ is not in W c for c 
strictly larger than 0, the proof is completed. □ 

The following is an outline for the rest of the article. We show that all solutions 
to ([7]) have bounded first and second spatial derivatives. This implies that all of 
their first partial derivatives are bounded (the time derivative is controlled by 1(2)1). 
Using the fact that ^ is autonomous in time, time translations of solutions are also 
solutions. We therefore construct a sequence of solutions {uu} to Cauchy problems 
started at t = 0, T\, T%, ... which tend to / + as t — ► +co, but their initial conditions 
tend to /_ as k — * oo. By Ascoli's theorem, this sequence converges uniformly on 
compact subsets to a continuous eternal solution. 
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3. Integral equation formulation 

In order to estimate the derivatives of a solution to |7jl . it is more convenient to 
work with an integral equation formulation of ([7]). This is obtained in the usual 
way. 



d_ 

Of 



du 
~di 

A ] u 



Au- u 2 + <fi 

-u 2 + (j) 



— - A 

dt 



■-U 2 ) 



(8) 

u(t, x) 



H(t,x-y)U(y)dy 



J-oo 



H(t - s,x-y) ((f>(y) - u 2 (s,y)) dyds, 



where H(t,x) — e~~ is the usual heat kernel. 

Calculation 4. We begin by estimating the first derivative of u for a short time. 
Let T > be given, and consider < t < T. The key fact is that J H(t, x)dx — 1 
for all t. Using © 



dx 



< 



< 



< 



< 



dU 




dx 


+ / 


dU 


+ /' 


dx 


oo -'0 


dU 




dx 


oo "'O 


dU 
dx 


+ T 

OO 



H(t - s, x - y) (<j){y) - u 2 (s, y)) dy ds 



./-oo 



d_ 

dx 

d_ 

dy 

. I ' dd> du 
H{t -^- y) {-d-y~ 2u dy- 



(H(t - s,x - y)) (<j>{y) - u 2 {s,y)) 



dy ds 



dy ds 



dx 



2\\u\ 



/' 


du 


/o 


dx 



ds. 



This integral equation fence is easily solved to give 



du 




dU 


+ T 







<- ( 






CO 


dx 


OO \ 


dx 


OO 


dx 





2tmax{j|/ + 



I/-IU} 



K,e K * T . 



Calculation 5. With the same choice of T as above, we find a bound for the 
second derivative in the same way: 



d 2 u 



dx 2 



< 



< 



d 2 U 


+ T 

CO 


d 2 4> 




dx 2 


dx 2 


+ / 

oo JO 


d 2 U 


+ T 

CO 


d 2 $ 




dx 2 


dx 2 


+ / 

oo JO 



d ( du 
dy V dy 



ds 



di 



dx 



+ 2\\u\\ 



dx 2 



ds 



for some Kz which depends on U, <fi, and T. 
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Calculation 6. Now, we extend Calculation [5] to handle t > T, 



d 2 u 


< 

oo 


d 2 U 




d 2 r T r°° 


dx 2 


dx 2 


+ 

OO 


dx 2 J J_ oc 



H(t -s,x- y)(<(>(y) - u 2 (s,y))dyds 



dx 2 

< K 3 e K2T + K 4 

< K 3 e K2T + K 5 

< K 3 e K * T + K & . 



H(t - s,x- y){<j>(y) - u 2 (s,y))dyds 



dx 2 



H(t - s,x) 



(||0||i + ||« 2 ||i)ds 



1 



■As + K' 



1 1 



i 



zds 



T s V s 

1 



tVt tVt 



hence there is a uniform upper bound on 
conditions, d>, and T. 



d- „ 



().r- 



which depends only on the initial 



Lemma 7. Let f 6 C 2 (K) be a bounded function with a bounded second derivative. 
Then the first derivative of f is also bounded, and the bound depends only on ||/||oo 
and ||/"||oo- 



Proof. The proof is elementary. The key fact is that at its maxima and minima, / 
has a horizontal tangent. From a horizontal tangent, the quickest /' can grow is at 
a rate of ||/"||oo- However, since / is bounded, there is a maximum amount that 
this growth of /' can accrue. Indeed, a sharp estimate is 



l/'lloo < v/Sll/lloolLfllc 



□ 



Using the fact that u is bounded, Lemma[7]implics that the first spatial derivative 
of u is bounded. By ([2]), it is clear that the first time derivative of u is also bounded. 

Lemma 8. As an immediate consequence of Lemmas&and^ the action integral 

+ -u 3 (t, x) — u(t, x)<p(x)dx 
o 

is bounded. Therefore, the solutions to the Cauchy problem (|7J| all tend to limits as 
t — > oo (Corollary 7 of ). By Lemma\^ we conclude that they all tend to the 
common limit of /+ when c > 0. 

Proof. The latter two terms are bounded due to the fact that u lies in L 1 D L°°(M) 
for all t. The bound on the first term comes from combining the fact that u and 
its first two spatial derivatives are bounded with the asymptotic decay of f± , and 
is otherwise straightforward (use L'Hopital's rule). □ 



A(u(t)) = 



,oo 1 


du 


Loo 2 


dx 
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4. Construction of an eternal solution 

Let 

U k (x) = (1 - 2- k - 1 )g 1/k (x) + 2- k - 1 f + (x), for k > 

noting that Uk — > f— as k — > oo. Since U k is a convex combination of /+ and 
gi/k, it follows that Uk € Wi/fc for all k. Also, since /+ and /_ have bounded first 
and second derivatives, the {Uk} have a common bound for their first and second 
derivatives. 

Now consider solutions to the following set of Cauchy problems 



(9) 



Uk(T k ,x) = U k (x). 



We choose T k so that for all k > 0, u k (0, 0) = u (0, 0). We can do this using the 
continuity of the solution and Lemma [8] As k — * oo, solutions are started nearer 
and nearer to the equilibrium /_, so we are forced to choose Tk — * — oo as k — > oo. 

It's clear that each solution Uk is defined for only t > Tk- However, for each 
compact set Scl 2 , there are infinitely many elements of {uk} which are defined on 
it. The results of the previous section imply that {uk} is a bounded, equicontinous 
family. As a result, Ascoli's theorem implies that {u k } converges uniformly on 
compact subsets to a continous u, which is an eternal solution to @- 

Our constructed eternal solution will have the value u(0,0) = Uq(0,0), which is 
strictly between / + and /_. As a result, the eternal solution we have constructed is 
not an equilibrium solution. By Lemma [SJ it is a finite energy solution, so it must 
be a heteroclinic orbit connecting /_ to /+. 
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